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Abstract 

Interfaces play an essential role in phonon-mediated heat conduction in solids, impacting appli¬ 
cations ranging from thermoelectric waste heat recovery to heat dissipation in electronics. From 
a microscopic perspective, interfacial phonon transport is described by transmission and reflection 
coefficients, analogous to the well-known Fresnel coefficients for light. However, these coefficients 
have never been directly measured, and thermal transport processes at interfaces remain poorly 
understood despite considerable effort. Here, we report the first measurements of the Fresnel trans¬ 
mission coefficients for thermal phonons at a metal-semiconductor interface using ab-initio phonon 
transport modeling and a thermal characterization technique, time-domain thermoreflectance. Our 
measurements show that interfaces act as thermal phonon filters that transmit primarily low fre¬ 
quency phonons, leading to these phonons being the dominant energy carriers across the interface 
despite the larger density of states of high frequency phonons. Our work realizes the long-standing 
goal of directly measuring thermal phonon transmission coefficients and demonstrates a general 
route to study microscopic processes governing interfacial heat conduction. 


^ These authors equally contributed to this work. 
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Transport across interfaces in heterogeneous media is a fundamental physical process 
that forms the basis for numerous widely used technologies. For example, the reflection 
and transmission of light at interfaces, as described by the Fresnel equations, enables wave- 
guiding with hber-optics and anti-reflection coatings, among many other functionalities. 
Interfaces also play an essential role in phonon-mediated heat conduction in solids.^ Material 
discontinuities lead to thermal phonon reflections that are manifested on a macroscopic scale 
as a thermal boundary resistance (TBR), also called Kapitza resistance, Rk, that relates the 
temperature drop at the interface to the heat flux flowing across it. TBR exists at the 
interface between any dissimilar materials due to differences in phonon states on each side 
of the interface.^ Typical interfaces often possess defects or roughness which can lead to 
additional phonon reflections and hence higher TBR. 

TBR plays an increasingly important role in devices, particularly as device sizes decrease 
below the intrinsic mean free paths (MFPs) of thermal phonons.^ At sufficiently small length 
scales, TBR can dominate the total thermal resistance. For instance, the effective thermal 
conductivity of a superlattice can be orders of magnitude smaller than that of the constituent 
materials due to high TBR.®^ This physical effect has been used to realize thermoelectrics 
with high efficiencjJ® and dense solids with exceptionally low thermal conductivity.^ On the 
other hand, TBR can lead to signihcant thermal management problemd^^^^ in applications 
such as and high power electronics.^^^^ 

Thus both scientihcally and for applications, a fundamental understanding of thermal 
transport across solid interfaces is essential. In principle, Fresnel transmission coefficients 
can also be used to provide a microscopic description of thermal phonon transport at in¬ 
terfaces owing to the similarities between photons and phonons. However, despite decades 
of work, the microscopic perspective of heat transport across interfaces remains poorly de¬ 
veloped compared to that available for photons. Today, interfaces are most often studied 
using macroscopic measurements of TBR or thermal conductivity. For example, numerous 
works have studied interfacial thermal transport by observing the temperature dependence 
of the thermal conductivitj^ or interface conductance, G = qj^ ] 3 y correlating 

changes in bonding strength and interface conduct ance.*^^*^ However, these experimental 
approaches provide limited information about the transmission coefficients because the ob¬ 
servable quantities are averaged over all phonons and thus obscure the microscopic processes 
of interest. 
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As a result, considerable uncertainty exists as to the values of the phonon transmission 
coefficients at solid interfaces. Commonly used analytical models for the transmission co¬ 
efficients include the gray model, in which the transmission coefficients are a constant for 
all phonon modes,l^®^ the acoustic mismatch modeP^® and the diffuse mismatch model 
(DMM),®® yielding incompatible predictions for the transmission coefficients. Further, 
their predictions are in poor agreement with experiment, and none are able to account 
variations in the atomic structure of actual interfaces. Atomistic methods such as molec¬ 
ular dynamicd^^®^ and atomistic Green’s functiond^i®^ have been extensively applied to 
obtain the transmission coefficients, predicting the phonon transmission coefficients to de¬ 
crease with increasing phonon frequency. However, no direct experimental verihcations of 
these predictions have been reported. Prior experimental observations have reached con¬ 
flicting conclusions, with one reporting that the transmission coefficients follow the trend of 
atomistic calculation^^ while another reaching the opposite conclusion that high frequency 
phonons are the dominant heat carriers across the interface.®^ Therefore, despite consid¬ 
erable experimental and theoretical study, the phonon transmission coefficients at actual 
interfaces remain unclear. 

Here, we report the hrst measurements of the Fresnel transmission coefficients for thermal 
phonons at a solid interface. Our approach is based on applying our recent advances in ab- 
initio phonon transport modeling based on the phonon Boltzmann transport equation (BTE) 
to interpret measurements from the time-domain thermoreflectance (TDTR), a widely-used 
optical experiment in the thermal sciences. Our measurements reveal that low frequency 
phonons are nearly completely transmitted across solid boundaries while high frequency 
phonons are nearly completely reflected, leading to an interfacial heat flux distribution 
that is dominated by low to mid frequency phonons. Further, our approach demonstrates a 
general route to directly experimentally study the microscopic transport processes governing 
interfacial heat conduction. 


TDTR WITH AB-INITIO PHONON TRANSPORT MODELING 

Our measurement is based on the TDTR experiment, an optical pump-probe technique 
that is used to characterize thermal properties on micron length scales. In this experiment, 
a femtosecond pulsed laser beam is split into a pump and a probe beam. The pump pulse 
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train is modulated at a frequency from 1 to 15 MHz to enable lock-in detection, and is then 
used to impulsively heat a metal him coated on the sample. The transient temperature 
decay at the surface is detected as a change in optical rehectance by the probe beam.H3 
Typically, this temperature decay curve is htted to a standard multilayer heat diffusion 
model based on Fourier’s law with the substrate thermal conductivity and metal-substrate 
interface conductance as htting parameters.^BSI This method is now a standard thermal 
metrology technique. 

Recently, considerable interest has focused on quasiballistic heat conduction in TDTR 
in which a lack of phonon scattering in the substrate leads to the breakdown of Fourier’s 
law. Many experimental reports have demonstrated clear evidence of quasiballistic heat 
transport in different material systems in the form of thermal properties that appear to 
deviate from their bulk values.®® In this work, we interpret these effects as fundamentally 
originating from the non-equilibrium phonon distribution emerging from the interface. As 
illustrated in Fig. [^a), when MFPs are much shorter than the characteristic length scale 
of the thermal gradient, information about the phonon distribution at the interface is lost 
due to scattering. On the other hand, if some phonons have sufficiently long MFPs, the 
non-equilibrium distribution penetrates into the substrate and affects the resulting heat 
conduction, thereby providing direct information about the spectral phonon distribution at 
the interface. 

While the required data is straightforward to obtain, the key to extracting the trans¬ 
mission coefficients is to rigorously interpret the data with the ab-inito phonon transport 
modeling based on the BTE. It is this step that has long impeded efforts to study interfaces 
due to the extreme cost of solving the BTE for the TDTR experiment. A number of sim¬ 
plified model^®® have been proposed to explain these observations. However, all of these 
models make various approximations that limit the information that can be extracted from 
the measurements. 

In this work, we overcome this challenge using two recent advances we reported for 
rigorously solving the spectral BTE under the relaxation time approximation (RTA), with 
no additional simplifications, that yield a factor of 10^ speedup compared to existing methods 
and allows the first ab-initio phonon transport modeling of TDTR free of artificial parameters 
or simplifications of the phonon dispersion. First, we have obtained an analytical solution of 
the spectral BTE in a semi-infinite substrate subject to an arbitrary heating profile.® Second, 
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we have employed a series expansion method to efficiently solve the spectral, one-dimensional 
(ID) BTE in a hnite layer, suitable for the transducer hlm.l^In this work, these two solutions 
are combined using a spectral interface conditioii^ that expresses the conservation of heat 
flux at each phonon frequency. Following the conclusions of prior computational work^^^*^ 
that the energy transmission at the interfaces considered here is elastic, we enforce that 
phonons maintain their frequencies as they transmit through the interface; direct electron- 
phonon coupling and inelastic scattering is neglected (see Supplementary Information). 

The details of the calculation are in the Supplementary Information. The only inputs to 
our calculation are the phonon dispersions and lifetimes, calculated using density functional 
theory (DFT) with no adjustable parameters by Jesiis Carrete and Natalio Mingo, and the 
only unknown parameter is the spectral transmission coefficients across the interface. We 
adjust the spectral transmission coefficients to obtain the best £t of the simulated surface 
temperature decay to the experimental measurement. The reflection coefficients are specihed 
by energy conservation once the transmission coefficients are known.^^ 


MEASUREMENTS OF PHONON TRANSMISSION COEFFICIENTS 

We demonstrate our transmission coefficient measurements by performing a TDTR mea¬ 
surement of an A1 him on Si substrate with the native oxide removed by Hydrohuoric acid 
prior to A1 deposition, yielding a clean interface. The TFM image in Fig. [^b) shows the 
interface thickness is less than 0.5 nm. The amplitude and phase of a typical signal from the 
lock-in ampliher are given in Figs, [^c) & (d). To begin, we follow the typical procedure of 
htting data with the standard heat diffusion model to extract thermal conductivity k and 
interface conductance Uj— We obtain G ~ 280 MW/m^-K and k ~ 140 W/m-K, which 
is in good agreement with prior works and literature values for the thermal conductivity of 
Si ]46l5^ 

Due to the good agreement, this measurement is typically taken as evidence that phonon 
MFPs in silicon are small compared to the thermal length ~ 1 pm, where a is the 
thermal diffusivity of silicon and /o is the modulation frequency. However, a number of 
ab-initio calculations, by different groups and with different computational packages, clearly 
show that a substantial amount of heat is carried by phonons with MFPs exceeding 1 
pm.l ^^ | 68 l This prediction has recently been experimentally confirmed by Cuffe et al using 
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FIG. 1. Measurements and simulations of TDTR experiments on Al/Si with a clean 
interface, (a) Schematic of the principle underlying the measurement of transmission coefficients. 
If the characteristic length scale of the thermal transport is much longer than the phonon MFPs, 
information about the interfacial distribution is lost due to strong scattering. If some MFPs are 
comparable to the thermal length scale, the non-equilibrium distribution at the interface propa¬ 
gates into the substrate where it can be detected, (b) TEM image showing the clean interface 
of a Al/Si sample with the native oxide removed. The interface thickness is less than 0.5 nm. 
Experimental TDTR data (symbols) on this sample at T = 300 K for a modulation frequency 
/ = 10.2 MHz along with the (c) amplitude and (d) phase fit to the data from the BTE simula¬ 
tions (dashed lines), demonstrating excellent agreement between simulation and experiment, (e) 
Apparent thermal conductivity and (f) apparent interface conductance of the experiments (sym¬ 
bols) and BTE simulations (lines) versus modulation frequency. These quantities are extracted by 
fitting the data and simulations to a thermal diffusion model. The pairs of solid lines denoted BTE 
bounds correspond to the uncertainty in the measured transmission coefficients plotted in Eig. 
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thermal measurements on variable thickness silicon membranes.^ This fact implies that 
quasiballistic transport should be readily observable in a typical TDTR experiment on Si, 
despite the seemingly correct thermal properties measured. 

This apparent contradiction can be resolved by considering the spectral prohle of the 
transmission coefficients. Our BTE calculations reveal that the measurements from TDTR 
on this sample strongly depend on the spectral prohle of the transmitted phonon spectrum, 
a dependence that does not occur in the heat diffusion regime. This dependence allows the 
transmission coefficients to be directly obtained from the TDTR data by adjusting them to 
ht the surface temperature data from the lock-in ampliher in a procedure exactly analogous 
to that used to measure bulk thermal properties. 

We performed the htting of transmission coefficients by adjusting them until the simulated 
surface temperature curve and experimental data from the lock-in ampliher matched each 
other. This comparison for all samples can be found in Supplementary Information. To 
more compactly report the data, we further process the BTE results by htting the simulated 
surface temperature decay curves to the same heat dihusion model used in the experiments to 
extract the apparent thermal conductivity and interface conductance at diherent modulation 
frequencies. If these two parameters at each modulation frequency match, then the htting 
curves will also match, enabling a compact comparison of the simulation and experimental 
data sets. However, we emphasize that the use of the heat dihusion model is for ease of 
comparison only and was not used in the transmission coefficient measurement. 

An example of the measurement process for the data shown in Figs, [^e) & (f) is given 
in Fig. Prior work^^^^used a constant transmission coefficient prohle that explained the 
apparent interface conductance. However, we hnd that this prohle predicts a modulation- 
frequency dependent apparent thermal conductivity that becomes as low as 100 W/m-K, 
in strong disagreement with experiment. Similarly, other commonly used models such as 
the DMM predict the wrong trend of thermal conductivity and interface conductance as a 
function of modulation frequency. (See Supplementary Information) 

The only transmission coefficient prohle that is able to simultaneously match the appar¬ 
ent thermal conductivity and interface conductance in Figs, [^d) & (e) from the experiment 
is shown in Fig. [^a). The hgure shows that the transmission coefficient from Si to A1 for 
longitudinal phonons, Tsi-s.Ai(i^), starts at unity, its maximum possible value and decreases 
steadily to near zero for high phonon frequencies. The transmission coefficient prohles for 
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FIG. 2. Transmission coefficients and spectral heat finx at the interface. Transmission 
coefficients of longitudinal phonons Tgi_^Ai(w) (blue shaded region) versus (a) phonon frequency 
and (b) phonon wavelength, along with the constant transmission coefficient profile (green dashed 
line) that gives the same interface conductance as the measured Tsi- 5 .Ai(w). The boundaries of the 
shaded region correspond to the BTE bounds in Figs, [^e) and (f). Our measurement shows that 
low frequency phonons have a much higher transmission probability across the interface than high 
frequency phonons, (c) Spectral heat flux with the measured (blue shaded region) and constant 
(green dashed line) transmission coefficient profiles across the interface versus phonon frequency, 
(d) Normalized accumulative interface conductance with the measured (blue shaded region) and 
constant (green dashed line) transmission coefficient profile versus phonon wavelength. Contrary 
to the prediction of the constant transmission coefficient profile, low-frequency/long-wavelength 
phonons carry a significant amount of heat across the interface. 

the other polarizations have similar shapes, and so thronghont the paper we plot only the 
longitndinal transmission coefficients for simplicity. The transmission coefficients from A1 to 












Si, TAi-!.si(i^) are calculated by satisfying the principle of detailed balance; the relationship 
between Tsi_s.Ai(a;) and TAi-).si(a;) reflect the differences in density of states and group ve¬ 
locity between the two materials. The transmission coefficients for each side of the interface 
and for the other polarizations are given in the Supplementary Information. 

Our measured transmission coefficient prohle thus indicates that low frequency, long 
wavelength phonons are transmitted to the maximum extent allowed by the principle of 
detailed balance, while high frequency, short wavelength phonons are nearly completely re¬ 
flected at the interface. This result is consistent with many prior works in the literature. 
First, our result qualitatively reproduces the physically expected behavior that transmission 
coefficient increases with decreasing phonon frequency because phonons with sufficiently 
long wavelength do not perceive the atomistic disorder at the interface. The absence of 
acoustic reflections in our TDTR signal conhrm the good acoustic match between the two 
materials and support our hnding that long-wavelength phonons are mostly transmitted at 
the interface. Second, the measured transmission coefficient prohle agrees with the experi¬ 
mental studies of polycrystalline silicon by Wang et al^ which suggested that transmission 
coefficient should decrease with increasing frequency. Finally, the measurement also agrees 
with numerous molecular dynamics and atomic Green’s function calculations, essentially all 
of which predict a decreasing transmission coefficient with increasing phonon frequency.^^^®^ 
Our work is thus able to provide direct experimental conhrmation of these predictions for the 
hrst time while eliminating other possibilities for the transmission coefficients that appear 
in the literature. 

Using this transmission coefficient prohle, we plot the spectral interfacial heat hux versus 
phonon frequency and accumulative heat hux versus phonon wavelength in Figs. ic) & (d). 
respectively. In contrast to the prediction by the constant transmission coefficient prohle, 
our results show that most of interfacial heat hux is carried by low to mid frequency phonons, 
with the contribution from high frequencies strongly reduced due to their small transmission 
coefficients. This observation highlights the importance of our ab-initio phonon transport 
modeling approach, as prior works neglected the contribution of these low-frequency phonons 
to interfacial heat hux.^ In fact, we are unable to explain the magnitude of the observed 
interface conductance without the contribution of phonons of frequency less than 3 THz. 
Similarly, we hnd that we can only explain the measurements using the exact phonon dis¬ 
persion for A1 computed from DFT; simple dispersion relations such as Debye model cannot 
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FIG. 3. TDTR measurements on Al/Si at different temperatures and Al/SiGe. (a) Ap¬ 
parent thermal conductivity and (b) apparent interface conductance versus modulation frequency 
from experiments (symbols) and simulations (shaded regions) for A1 on Si with a clean interface 
at 300 K, 350 K and 400 K. (c) Apparent thermal conductivity and (d) apparent interface con¬ 
ductance versus modulation frequency from experiments (circles: this work; triangle: Ref. 50) and 
simulations (shaded regions) for A1 on SiGe with a clean interface. The magnitude and trend of 
the experimental data are reproduced using the same transmission coefficient profile as in Fig. 

explain the data because they underestimate the contribution of low frequency phonons to 
thermal transport. 

Since the energy transmission at the interfaces is considered elastic, the transmission 
coefficient prohle in theory should be independent of temperature. To conhrm the robustness 
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of the measured transmission coefficients, we conducted experiments at 350 K and 400 K 
and compared the experimental results with the calculations using the same transmission 
coefficient prohle measnred at 300 K. As shown in Figs. (a) & (b), the transmission 
coefficient prohle measured at 300 K yields excellent agreement between simulations and 
experiments at higher temperatures. 

To further support our measurements, we additionally measure the transmission coeffi¬ 
cients for A1 on SiGe. While this material has an additional point defect scattering mecha¬ 
nism compared to pure Si, we expect the transmission coefficients to be nearly the same given 
that the host lattice is nnchanged. Figs, [^c) & (d) plot the apparent thermal condnctivity 
and interface conductance, demonstrating that our BTE modeling with the exact transmis¬ 
sion coefficient prohle shown in Fig. |^a) agrees well with the measured apparent thermal 
condnctivity and interface conductance for Al/SiGe. This result conhrms that the measured 
transmission coefficients for Si and SiGe snbstrates are indeed the same. Importantly, the 
apparent thermal conductivity of SiGe dihers from its actual thermal conductivity, mea¬ 
sured by Thermtest as 50.7 ± 0.5 W/m-K using the transient plane sonrce method on a 
bnlk sample. This discrepancy highlights the strong inhnence of the transmission coefficient 
prohle on the measnrements from TDTR. 

So far, the measurements indicate that low freqnency phonons are nearly completely 
transmitted throngh a clean interface. We next seek to determine what types of interfacial 
disorder can rehect these modes. We condncted additional measurements for A1 on Si with 
a native oxide layer (thickness ~ 1 nm as shown in a TEM image in Fig. |^a)) and Si with 
thermally grown oxide layer (thickness ~ 3.5 nm as shown in a TEM image in Fig. |^b)). 
Since the oxide layers are sufficiently thin to neglect their thermal capacitance, we can treat 
them as part of the interface that modihes the transmission coefficient prohle in onr current 
BTE model. 

By again htting the BTE results to the measurement as in Figs. |^c) & (d), we are 
able to hnd the transmission coefficient prohles for these two cases as shown in Figs. |^e) 
& (f). Gompared to a clean interface, the transmission coefficients for A1 on Si with a 
native oxide are rednced for all phonon modes. However, low-freqnency phonons experience 
a larger reduction in transmission coefficient than do high-frequency phonons, which have 
transmission coefficients close to zero even at a clean interface. When the ronghness of the 
interface increases with a thicker oxide layer, the transmission coefficient keeps decreasing 
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FIG. 4. Transmission coefficient measurements of Al/Si with different types of inter¬ 
faces. TEM images showing the Al/Si sample with (a) native oxide layer (thickness ~ 1 nm) 
and (b) thermally grown oxide layer (thickness ~ 3.5 nm). (c) Apparent thermal conductivity 
versus modulation frequency and (d) apparent interface conductance versus modulation frequency 
of experiments (symbols) and simulations (shaded regions) for A1 on Si with native oxide layer and 
Si with thermally grown oxide layer. The lines give the upper and lower bounds of the BTE simu¬ 
lations used to determine the transmission coefficients. The corresponding transmission coefficient 
profiles versus (e) phonon frequency and (f) phonon wavelength show that as the interface gets 
rougher, low frequency/long wavelength phonons are more likely to be reflected. 
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and more phonons, especially those with wavelengths between 1 and 2 nm, are reflected at 
the interface. Therefore, our measurements show that phonons with wavelength shorter than 
the interface roughness are more likely to be reflected by the interface than long-wavelength 
phonons, and as the interface gets rougher, a larger fraction of the phonon spectrum is 
affected by the interface. 


DISCUSSION 

Our work has considerable implications for thermal metrology and technological applica¬ 
tions. First, TDTR is a widely used characterization tool to measure the thermal properties 
of thin Aims and bulk materials. However, this work demonstrates that in the quasiballistic 
regime, the apparent thermal conductivity and interface conductance measured by TDTR 
may not reflect bulk properties because both are strongly affected by the transmission coef¬ 
ficient profile. The apparently correct measurement of silicon thermal conductivity appears 
to be a fortuitous cancellation of two factors: the high transmission coefficient of low fre¬ 
quency phonons leads to an increased contribution to heat flux that offsets the deviation 
from Fourier’s law that occurs due to a lack of scattering. If these two factors were not 
balanced, the apparent thermal conductivity of Si would not coincide with the bulk value. 
In fact, these factors are not balanced in SiGe, resulting in an apparent thermal conductiv¬ 
ity that is substantially smaller than the actual value. Additionally, the apparent thermal 
interface conductance measured in SiGe deviates from the value measured in Si despite the 
fact that the transmission coefficients are exactly the same. These observations are further 
evidence that the apparent thermal properties extracted in a conventional TDTR inter¬ 
pretation approach do not necessarily reflect the actual physical properties of the materials. 
Although this situation is undesirable from a thermal metrology perspective, our work shows 
that TDTR is capable of providing considerable microscopic detail about thermal phonons 
provided the measurements can be properly interpreted. 

Second, our measurements show that the spectral profile of transmission coefficients is 
essential to understanding thermal transport across interfaces. Due to a lack of knowledge 
about interfaces, the phonon transmission coefficients are often taken to be a constant or 
calculated using the DMM. However, this work shows that both of these simple models 
are incapable of explaining the experimental measurements. Although the gray model and 
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measured transmission coefficient profiles yield the same apparent interface conductance, 
the spectrum of heat transmitted across the interface is completely different, resulting in 
opposite conclusions regarding the primary heat carriers across the interface. Therefore, 
considering the spectral transmission coefficient prohle across an interface is essential to 
accurately describing thermal phonon transport. 

Third, our work provides strong evidence that elastic transmission of phonons across an 
interface is the dominant energy transmission mechanism for materials with similar phonon 
frequencies. Our BTE model does not incorporate electrons or inelastic scattering yet is able 
to explain all of the measurements we performed. We conclude that inelastic transmission 
and coupling between electrons in metals and phonons in semiconductors have little influence 
on the energy transport for the materials considered here. 

Fourth, our results demonstrate that disorder at interfaces plays an important role in 
the spectral content of the heat transmitted through the interface and provides strategies 
to alter interface conductance. For instance, in applications like LFDs where the heat dissi¬ 
pation rate across interfaces is to be enhanced, the key to increasing interface conductance 
is to minimize the reflection of high frequency phonons by reducing defects; low frequency 
phonons are likely to be mostly transmitted already. On the other hand, the strong frequency 
dependence of the transmission coefficients can be exploited to create thermal phonon £1- 
ters to selectively remove parts of phonon spectrum, analogous to optical long-pass filters. 
Phonons with wavelength much longer than the characteristic roughness of an interface are 
more likely transmitted through the interface while short-wavelength phonons are mostly 
reflected. The thermal phonon spectrum responsible for heat condnction can thus be ma¬ 
nipulated by controlling the atomistic ronghness of an interface. 

Finally, onr work exemplifies the powerfnl insights that can be obtained throngh ab- 
initio phonon transport modeling, free of artificial parameters or arbitrary divisions of the 
phonon spectrnm into diffnsive and ballistic modes. Most of the previons stndies focnsed on 
developing simplified thermal models to explain the modulation-dependent measnrements for 
alloys snch as SiGe bnt did not consider the contradictory observations in silicon. Obtaining 
the transmission coefficients from measnrements that appear to be in good agreement with 
diffnsion theory would be challenging withont the predictive calcnlations reported here. Onr 
work thus demonstrates a general approach that can be used to advance onr knowledge of 
thermal phonons to a level on par with that of photons. 
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METHODS 


Sample preparation Commercial high-purity natural Si (100) wafer and Si-Ge (1.5-2 
at % Ge) wafer (100) from MTI Gorp. were used in the experiments. Before coating A1 
on the samples, three different surface conditions of the samples were prepared. First, the 
native oxide was removed with buffered HF acid to obtain a clean surface of Si and SiGe. 
After etching, the samples were immediately put into a vacuum chamber for A1 deposition. 
Second, the native Si02 layer was left in place. No further treatment was taken for this 
condition before A1 deposition. Finally, a thermally grown Si02 layer as fabricated by 
putting the Si samples into a tube furnace for three hours. The thickness of the native Si02 
layer and thermally grown Si02 layer was measured by ellipsometry and TEM to be ~ 1 nm 
and ~ 3.5 nm, respectively. A thin him of A1 was deposited on all samples using electron 
beam evaporator. The thickness of the A1 transducer layer was 70 nm, measnred by atomic 
force microscopy. 

TDTR measurements The measnrements are taken on two-tint TDTR. The details 
are available in Ref. 64. The probe diameter is 10 /rm and the pnmp diameter is 60 /xm. 
Both beam sizes are measnred using a home-bnilt two-axis knife-edge beam prohler. With 
60 /im pump heating size, the heat transfer problem can be treated as one-dimensional. All 
the measurements at T = 300 K are performed nnder ambient conditions, and the additional 
measnrements at T = 350 and 400 K are performed in an optical cryostat (JANIS ST-100) 
under high vacnnm of 10“® torr. 

TEM images The TEM samples were prepared by standard FIB lift-out technique in 
the dual beam FE-SEM/FIB (FEI Nova 600). To protect the top surface, a Pt layer with 
thickness ~ 300 nm was deposited with electron beam evaporation followed by another Pt 
layer with thickness ~ 3-4 /rm by Ga ion beam. The lamella was cut parallel to the chip edge 
which was aligned to the wafer flat edge dnring initial cntting in TDTR sample preparation. 
As a result, the cutting snrface normal was along (110) direction and all the TEM images 
were taken parallel to the Si (110) crystallographic zone axis. High resolution transmission 
electron microscopy (HRTEM) analyses were carried ont in a FEI Tecnai TF-20 TEM/STEM 
at 200 kV. To avoid damage from the high energy electron beam, the beam exposure on 
region of interest was minimized especially at high magnihcation dnring operation. 

Ab-initio modeling Thermal transport in TDTR experiments, assuming only cross- 
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plane heat conduction, is described by the spectral, one-dimensional (ID) BTE under RTA.I^ 
Efficiently solving this equation for the conditions of the TDTR experiment required two of 
our recent computational advances. First, we have obtained an analytical solution of the 
spectral BTE in a semi-inhnite substrate subject to an arbitrary heating prohle.^ Second, 
we have employed a series expansion method to efficiently solve the spectral ID BTE in a 
hnite layer, suitable for the transducer hlm.l^Sl In this work, these two solutions are combined 
using a spectral interface conditiorP^ that expresses the conservation of heat flux at each 
phonon frequency. Following the conclusions of prior computational workd^^^^ that the 
energy transmission at the interfaces considered here is elastic, we enforce that phonons 
maintain their frequency as they transmit through the interface; direct electron-phonon 
coupling and inelastic scattering are neglected. The details of the calculation are in the 
Supplementary Information. The only inputs to our calculation are the phonon dispersions 
and lifetimes, calculated using DFT with no adjustable parameters by Jesiis Carrete and 
Natalio Mingo, and the only unknown parameters are the spectral transmission coefficients. 
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